Transport and localization of waves in one-dimensional disordered media : Random 

phase approximation and beyond 
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We report a systematic and detailed numerical study of 
statistics of the reflection coefficient (\R(L)\ 2 ) and its as- 
sociated phase (ff) for a plane wave reflected from a one- 
dimensional (ID) disordered medium beyond the random 
phase approximation (RPA) for Gaussian white-noise disor- 
der. We solve numerically the full Fokker-Planck (FP) equa- 
tion for the probability distribution in the (\R(L)\ 2 , 6(L))- 
space for different lengths of the sample with different "dis- 
order strengths". The statistical electronic transport prop- 
erties of ID disordered conductors are calculated using the 
Landauer four-probe resistance formula and the FP equation. 
This constitutes a complete solution for the reflection statis- 
tics and many aspects of electron transport in a ID Gaussian 
white-noise potential. Our calculation shows the contribu- 
tion of the phase distribution to the different averages and its 
effects on the one-parameter scaling theory of localization. 



PACS numbers: 72.15.-v, 72.15.Rn, 73.23.-b, 05.40. +j 



I. INTRODUCTION 

Wave propagation through ID random media have 
been studied for decades. Typical physical examples of 
these problems are electron transport in disordered con- 
ductors, light transport in random dielectric media, etc.. 
All previous studies agree that the transmission through 
a ID disordered medium decays exponentially with the 
length of the sample. For electronic systems, the conduc- 
tance decays exponentially with the length of the sample 
and the resistance increases exponentially with the sam- 
ple length. The main characteristic length scale in the 
problem is the localization length. Results of different 
phases of research in this field have been reviewed re- 
cently in several review articles |l]-|7j . 

Disordered ID systems are quite generic as regards 
their transport properties — disordered ID conductors 
are all alike, but every ordered ID conductor is ordered 
in its own way! Thus, Mott et. ai. |Q first showed that 
all the eigenstates of I D disordered quantum systems are 
exponentially localized. There are no truly ID disordered 
metals. For electronic systems quantum interference ef- 
fects are important at low temperatures when decoher- 
ence due to the inelastic process may be neglected. At 
low enough temperatures with static disorder, the trans- 
port properties of the electronic systems are highly fluc- 



tuating (i.e. sample specific) due to the different con- 
ductors with the same impurity concentration but with 
different microscopic arrangements of the impurities, can 
differ substantially in their transport properties. Sev- 
eral studies show that the root mean-squared-fluctuation 
is more than the average for length scales larger than 
the localization length, while for the length scales within 
the localization length (i.e. in the good metallic regime 
for quasi fD and higher dimensional systems) the fluc- 
tuations are finite and universal, the so called univer- 
sal conductance fluctuation (UCF) 0. These fluctua- 
tions make the resistance and the conductance non-self- 
averaging quantities. The resistance at low temperatures 
when quantum effects are important is non-Ohmic, i.e., 
non additive in series. The non-additive nature arises 
because of the non-local effects of the quantum wave am- 
plitudes associated with the electron in the conductor. 
One has to take proper care of the phases to add the two 
quantum resistances together. 

To calculate any meaningful quantity for a non-self- 
averaging quantity, like the resistance and the conduc- 
tance, one should know the full probability distribution of 
the same for every length of the sample. Several authors 
|^| |l2| have derived a Fokker-Planck equation for the full 
probability distribution of the resistance for the Gaus- 
sian white noise disorder. Exact results for the average 
resistance and the average conductance can be calculated 
analytically for the weak-disorder case using the Fokker- 
Planck equation. The approximation made in the weak- 
disorder case is that the phase of the complex reflection 
amplitude relative to the incident wave, or the relative 
phase of two non-ohmically additive quantum resistances, 
are distributed uniformly, which is the random phase ap- 
proximation (RPA). Then the problem can be solved ana- 
lytically to calculate the average resistance exactly for the 
Gaussian white noise potential. The analytical results 
show that the resistance and the conductance have log- 
normal distributions for large lengths. The average four- 
probe resistance increases exponentially with the length 
of the sample with a characteristic length scale called the 
localization length. The average four-probe conductance, 
however, diverges for all length scales. The divergence of 
the average conductance for all length scales is due to 
resonance (absolute transmission) states which dominate 
the probability distribution jl3],[l4| . We will discuss else- 
where the issues related to the divergence of the average 
conductance. 

It is difficult to get an average quantity by direct ana- 
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lytical calculation which includes the phase distribution, 
where the phase distribution depends on the strength of 
the disorder. However, the actual validity of the ran- 
dom phase approximation has not been studied in detail. 
There are issues like actual contribution of the phase dis- 
tribution to the different averaging processes that have 
not been systematically studied so far. 

We study numerically, the joint probability distribu- 
tion of the phase (9) and the reflection coefficient (r) 
of the complex amplitude reflection coefficient (R(L) = 
yfr exp(i9)) for a one-dimensional disordered conductor 
with the Gaussian white noise disorder, for different 
lengths of the sample and different strengths of the dis- 
order. This is mainly the solution of the Fokker-Planck 
equation, where the probability density is varying in the 
(r, 6>)-space and evolving with the sample length. Us- 
ing the "invariant imbedding" technique, a non-linear 
Langevin equation can be derived for the R{L). Then, 
using the stochastic Liouville equation for the probabil- 
ity evolution and the Novikov theorem to integrate out 
the stochastic aspect due to the random potential, one 
gets the Fokker-Planck (FP) equation in the (r, 0)-space 
for varying sample length L jf0|,[l2]]. Integrating 9 and 
r parts separately from the joint probability distribution 
P(r, 9), we have calculated the marginal probability dis- 
tributions for r and 9, respectively, for a given length 
and disorder strength of the sample. The Landauer for- 
mula |l5|,|l6| relates the resistance and the conductance 
directly to the reflection coefficient (r). Once the reflec- 
tion probability is known, the average electronic prop- 
erties of the system such as resistance, conductance etc. 
can be calculated easily using the Landauer formula. We 
have outlined the range of validity of the random phase 
approximation, for the parameters on which the phase 
distribution depends, namely the localization length and 
the wave vector of the incoming wave. Our calculation 
shows that £ > A = 2it /k is the limit where the random 
phase approximation is valid, where £ is the localization 
length and A is the wave length of the incoming wave (and 
for the electronic case we will consider A to be Fermi wave 
length) 

Knowing the probability distribution of the reflec- 
tion coefficient and using the Landauer four-probe resis- 
tance/conductance formula, we have calculated the prob- 
ability distribution for the resistance/conductance and 
find it to be log-normal, for all strengths of the white 
noise Gaussian disorder. We calculate the averages and 
the fluctuations of the reflection coefficient, the conduc- 
tance and the resistance and of logarithms of the con- 
ductance and the resistance. Lastly, we study the effects 
of the phase distribution on the one-parameter scaling 
theory of localization. We have shown that the one- 
parameter scaling theory holds quite well for large sample 
lengths, but the phase has a definite but small effect on 
it. To the best of our knowledge this is the first work 
where the joint probability distribution for the reflection 



coefficient and its phase for a ID Gaussian white-noise 
disordered conductor is calculated for all length scales 
and for different disorder strengths. The way we have cal- 
culated the different quantities here involves essentially 
no approximations. We will discuss works of others while 
discussing our results. We report here mainly the results 
of the electronic transport in disorder conductors, but 
some parts of the formalism apply equally well to the 
case of light transport in random dielectric media. 



II. METHOD OF CALCULATION 

Both the Schrodinger and the Maxwell equations can 
be transformed to the Helmholtz equation. Therefore, 
these two equations can be studied in the same frame 
work for the reflection and the transmission amplitude, 
since they transform to the same Helmholtz equation and 
have similar initial conditions. 

Consider the Schrodinger Equation, 



d 2 ib 



(1) 



where we have set h 2 /2m = 1, 
and the Maxwell equation, 



d 2 E u 2 r 

_ + _ [eo + e{x)]E =0 . 



(2) 



Transforming the Schrodinger and the Maxwell equation 
to the standard Helmholtz equation form, we obtain : 



d 2 u 
dx 2 



+ k 2 [I + n(x)] u =0, 



(3) 



where rj(x)=— V(x)/k 2 for Schrodinger wave, and rj(x)= 
and k 2 = ^-eo f° r the Maxwell wave, where eo is 
the constant dielectric background and e(x) is the spa- 
tially fluctuating part of the dielectric medium. 



A. The Langevin Equation for the complex 
amplitude reflection coefficient R(L) 

Consider a plane wave of wave vector k incident from 
the right side of the disordered sample of length L 
having the reflection amplitude R(L). The non- linear 
Langevin equation for the reflection amplitude for the 
plane wave scattering problem can be derived by the in- 
variant imbedding technique Jl2[ . The Langevin equation 
for the complex amplitude reflection coefficient is : 



dR(L) 
dL 



2ikR(L) 



R{L)f 



(4) 



with the initial condition R(L)=0 for L=0. 

The main idea is to use the Langevin equation to get 
a FP equation for the reflection probability density. For 
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the electronic system consider the reflection coefficient 
(or the transmission coefficient) of a free electron wave 
scattered by the random potential in order to calculate 
the resistance/ conductance of the system. The average 
properties of the resistance/conductance can be calcu- 
lated using the Landauer four-probe formula and the FP 
equation for the reflection coefficient. 



B. The Fokker-Planck equation 

To get the Fokker-Planck equation from the non-linear 
Langevin equation one has to get the probability 
density equation first and then the stochastic aspect due 
to the random potential has to be integrated out. The 
Langevin equation (EqJ^) can be solved analytically for 
the Gaussian white noise potential to get the FP equa- 
tion. The non-linear Langevin equation for R(L) in 
(Eq.^) is basically two coupled differential equations for 
the magnitude and the associated phase parts. 
Taking 



R(L) = v/r(Z)«Bp(i0(L)), 



(5) 



and substituting Eq.^| into Eq.^, and equating the real 
and the imaginary parts on both sides of the Eq^, one 
gets, 



dr 
dL 



kn(L)r^ (1 



(6) 



< n(L) >= and < T](L)ri(L') >= q5(L - L'). (10) 

Eq.(||) has terms like rjQ which are to be averaged out. 
For the Gaussian white noise disorder, the Novikov the- 
orem says : 



< 



SQ[y] 
6 V (L) 



> 



(11) 



After averaging out the disorder aspect in EqJ^ and writ- 



ing < Q(r, 
equation: 



> v = P(r 



one gets the Fokker-Planck 



dP(r, i 
df 



smtf— r i/ ^(l - r) + 



dr 



09 



V /2 



P{r,6) 



2fc£ 



dP(r, I 
W 



(12) 



Where I = L/£ and £ = (^qk 2 )^ 1 is the localization 
length. 

The Fokker-Planck equation, Eq.[l2] has all the infor- 
mation of the probability distribution of the reflection 
coefficient (r) and the associated phase (9) for different 
length scales of the sample and with varying disorder 
strengths. We will explore several aspects of the Eq|l2| in 
our study, which we are going to discuss in detail in later 
Sections. 



dL- 2k+ 2^ 



2 + cos 9{r 



1/2 



(7) 



Now, according to the van-Kampen lemma |17[ , these two 
stochastic coupled differential equations will produce a 
flow of the density Q(r, 9) in the (r, #)-space according to 
the stochastic Liouville equation with increasing length 
of the sample, i.e., Q(r, 9) is the solution of the stochastic 
Liouville equation: 



dQ{r,9,L) 

dL 



d ,^dr , d ,^d9 , 

where dr/dL and d9/dL are given by Eqsj6|and 0. Now, 
substituting the values of dr/dL and d9/dL in the Eq.|| 
one gets, 



dQ 
dL 



dr 



r x/2 (l-r) V (L)Q 



^ Vi{l)Q] 



~^-'l(r 1/2 +r-^)^[co S 9 V (L)Q}. (9) 

To get the Fokker-Planck equation, Eq.|| has to be av- 
eraged over the stochastic aspect, i.e., over all realiza- 
tions of the random potential. For the case of a Gaus- 
sian white noise potential, EqJ^ can be averaged out over 
the stochastic potential analytically using Novikov's Q 
theorem. For the Gaussian white noise potential: 



III. PARAMETERS OF THE PROBLEM 

The Fokker-Planck equation (Eqjl2[) has three param- 
eters to describe the problem fully: (1) the length of the 
sample L, (2) the localization length £, and (3) the inci- 
dent wave vector k. In the re-arranged form of the Eqjl2|, 
as it is written, it has effectively two parameters: I — L/£ 
and C = 2k£. Here I is a number which gives the length 
of the sample in units of the localization length and C 
is a number which fixes inverse of the disorder strength 
in terms of the wave vector of the incident wave and the 
localization length. Larger value of C implies that the lo- 
calization length is large, or the incoming electron energy 
is higher, or both, that is, the weak-disordered regime. 
Conversely, when C is small it means £ is small, or the 
incoming wave energy is small or both, that is, the strong 
disorder regime. 



IV. ANALYTICAL SOLUTIONS FOR THE 
RESISTANCE AND THE CONDUCTANCE IN 
THE RANDOM PHASE APPROXIMATION 
(RPA) 
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A. Resistance 

In the random phase approximation (RPA), which is 
valid for weak disorder and large incident electron en- 
ergies, one can write P(r,0) = (l/27r)P(r) i.e. P(r,0) 
factorizes, and 8 is uniformly distributed over 2ir. Con- 
sidering dP/dd = 0, the Fokker-Planck equation Eq.|l2] 
in r then becomes: 



dP{r) 
dl 



d_ 

Or 



r^-{l-r) 2 P{r) 



(13) 



( Here we have used the same symbol P(r) for the 
marginal probability density of r as for the joint density 
P(r,0) .) 

Now, from the Landauer four-probe resistance formula, 
the dimensionless resistance as a function of the reflection 
coefficient is: 



P(l) 



r(0 



1 



■(I)' 



(14) 



Making the transformation from P(r) to P(p) by us- 
ing P(p) — P(r)(dr/dp), one obtains the Fokker-Planck 
equation for the resistance: 



dP(p) d , d , 

-ir = d-p p{l+p) d- P p ^ 



(15) 



with the initial condition, P(p) = 8(p) for I = 0. 

Eq.|l5| has been derived earlier by several authors 
p9| , |l0yi2] . This equation has also been derived from the 
maximum entropy principle (MEP) by Mello and Kumar 
{§• 

The average of p n can be obtained without solving di- 
rectly Eq.|l5| as following. 
Let us define, 



P{p)p n dp 



(16) 



Multiplying EqjlJ by p n and integrating both sides of the 
equation for p from to oo, one gets a moment recursion 
equation for the average moments of the resistance, 



d Pn , . .2 

— - = n(n + l)p n + n p n -\- 



(17) 



Since the probability is always normalizable, the value of 
po will be: 



poo poo 
PO = / P(p)p° = / P(p) = 1 

J o Jo 



(18) 



Once we know the initial value po, then Eq.|l7] can be 
solved analytically for higher moments of the average re- 
sistance, 



Pi(0 = 5 ( e 

p 2 (0 = ^(2e 6i -6e 2i +4), 



PnQ) 



a n(n+l)l 



n > 1 



n(ri + l)!' 
and specifically, 

rms fluctuation of p ~ e 31 . 



(19) 
(20) 
(21) 

(22) 



The above expressions imply that the average resistance 
increases exponentially with the length of the sample and 
the rms fluctuation is greater than the average, indicat- 
ing that the resistance is not a self-averaging quantity. 
It is clear now why one has to consider the full prob- 
ability distribution to describe the disordered quantum 
resistors. Eq.^B] has also an analytical expression for the 
full distribution of p for the large I limit Jl(J, which is 
log-normal. 



B. Conductance 

If one solves the Fokker-Planck equation for < g n > 
(=< l/p n >), it has been shown by several workers and 
explicitly by Melnikov [[l9| that all the moments v of 
< g v > diverge for v > 1/2, for all length scales of the 
sample. The cause is the existence of the resonant states, 
i.e., states r ~ 0. Though, the probability is less for the 
resonant states, but it is finite. The existence of the 
resonance states have been discussed in detail in Ref. 

us 

V. AIM OF OUR WORK 

The main aim of our work is to evolve the full Fokker- 
Planck equation, (Eq.[l2|), numerically. Eq.|l2| has full 
informations of the the statistics of the reflection coeffi- 
cient and its phase for the case of Gaussian white noise 
disordered potential. We will study several aspects of 
the solution, and try to see the actual contribution of 
the phase, which goes beyond the random phase approx- 
imation. Important point to note here is that the nu- 
merical calculation part involves only the solution of the 
Fokker-Planck equation which has been derived analyti- 
cally without approximation. 



VI. NUMERICAL DETAILS 

We took r and 9 as Cartesian variables in a two- 
dimensional 50 x 50 grids. An explicit finite-difference 
scheme |2^j2^] was used to solve the FP equation. The 
von-Neumann stability criterion was checked and the 
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Courant condition for the used discrete iterative length 
was strictly maintained. A few results were also checked 
by using the rather time consuming implicit finite- 
difference scheme. 

A. Allowed error bars 

Error bars of the order of 10 -4 for r, 3 x 10~ 3 for 9, 
and 10 -12 for length I, were allowed for the whole range 
of numerical calculations. 

B. Initial probability distribution P(r,6) at L — 

The FP equation (|l|) poses an initial value problem. 
The initial probability distribution P(r, 0) at I = has 
to be specified, which will then evolve with the increase 
of the length of the sample. The Fokker-Planck equa- 
tion ( ]l2| ) is however singular at r=0. This causes a tech- 
nical problem for solving the equation numerically. To 
circumvent this problem, we have therefore taken an ini- 
tial (fixed) scatterer with r = .01 by putting a half-delta 
function potential peaked at I = which could be phys- 
ically understood as due to an initial impurity sitting 
at I = 0, or the contact resistance of the leads. By 
"fixed" we mean that it is fixed over all the realizations 
of the sample randomness. Phase distribution for such a 
weak delta-function potential will peak around +tt/2 or 
— it/2, depending on the sign of the delta-function poten- 
tial. Once ro = .01 is fixed, then the phase distribution 
has equal probability peak at +ir/2 and —tt/2. A fixed 
weak delta scatterer at the position I = will not change 
the gross statistics, except at very smaller length scales. 
We have kept this initial distribution same throughout 
the numerical calculations. We could not consider any 
smaller value of the initial- fixed- reflection coefficient (ro), 
or a lower cut-off to the r = singularity, for reason of 
convergence criterion of the numerical algorithm. An es- 
timate can be done for the initial cut-off length Iq of the 
sample (in terms of the localization length £) for this 
small r = ro. Taking analytical results for the weak dis- 
ordered case, one gets: 

Po = 7 -^ L ^ = l[exp(2l )-l} (23) 

then 

l = L /£ = 2ln[\±^] ~ hn(l + 2r ) ~ r = .01 (24) 
1 - ro 2 

This implies that the initial length is 1% of the localiza- 
tion length, throughout the numerical calculation. 

For numerical calculation, the delta- function has to be 
taken as the limit of a continuous function. In (r, &)- 
space for a physically reasonable initial probability dis- 
tribution, we have taken this as 



P(r, 9) l=0 = S(r - .01) [6(9 - tt/2) + 5(9 + tt/2)} , (25) 

where the delta functions are sharp Gaussians. 

Fig.|l| shows this initial probability distribution P(r, 9) 
nominally at I = 0, and the marginal distribution P(r) 
of the reflection coefficient and the marginal distribution 
P(9) of the phase 9. It should be emphasized again here 
that the initial probability distribution of the phase 9 
can be taken at any small enough length. However, the 
statistical properties of the system do not depend on the 
initial distribution except for very small lengths. 

C. Boundary conditions for r and 9 for any length L 

The unit step length of the discrete evolution is taken 
to be AZ = 10~ 6 . For every discrete evolution, we took 
P(r, 9) = for r > 1 and r < along r axis; and the 
boundary condition was taken as periodic along 9 axis, 
such that P{r,2ir + A9) = P(r,A9) for every discrete 
evolution. 

VII. RESULTS AND DISCUSSIONS 

A. Evolution of P(r, 0) with the length L for different 
disorder strengths 

We will consider the evolution of the full probabil- 
ity distribution P(r, 9) for different lengths I for the 
three main regimes of disorder strength — (1) weak, (2) 
medium, and (3) strong. 

1. P(r,0) for weak disorder : 2k£ = 100 (2k£ » I) 

Figs.^(a),(b),(c),(d) show typical distributions of the 
probability in (r, #)-space for different lengths I — 
1, 2, 5, and 10. Evolution shows clearly that the phase 
part is uniformly distributed. Evolution of the proba- 
bility is mainly along the amplitude (r) axis. Probabil- 
ity evolves with length and peaks near r = 1 for large 
lengths. These evolution pictures imply that the random 
(uniform) phase distribution holds for the weak disorder, 
or for the large values of the localization length. Later we 
show the nearly exact range for the disorder parameter 
(2fc£) where the random phase approximation is valid. 

2. P(r, 6) for medium disorder: 2k£ — 1 (2k£ ~ 1) 

Similarly, Figs.|| (a),(b),(c),(d) are as Fig.|| but for 
the disorder parameter 2fc£ — 1, i.e., medium disorder 
regime. These evolutions clearly show that the probabil- 
ity distribution P[r, 9) is quite complicated and has an 
asymmetric distribution in phase. The phase distribution 
peaks on one side (9 > tt) for large lengths. 
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3. P(r,6) for strong disorder : 2k£ = .001 (2fc£ << 1) 

Similarly, Figs.[|(a),(b),(c),(d) are as Fig.||, but for the 
disorder parameter 2fc£ = .001, i.e., the strong disor- 
der regime. Probability distribution is perfectly sym- 
metric centered at the phase 9 = tt. In the strong dis- 
order parameter regime the sample tries to behave as 
a perfect reflector and totally reflect back the incoming 
wave with opposite phase. Evolution pictures show that 
P(r, 9) peaks around r = 1 and is symmetrically around 
9 = tt for larger lengths. It will be shown later that the 
distribution does not change substantially with further 
increase of the disorder strength parameter, i.e. the dis- 
tribution is insensitive to the disorder parameter 2fc£ in 
that regime. 



B. Marginal distribution P(r) of the reflection 
coefficient (r) and Marginal distribution P{9) of the 
phase (6) for strong, medium and weak disorders 

Figs.|(a),(b),(c) show (i) the marginal probability dis- 
tribution P(r) when the phase((9) part of the P(r, 9) is 
integrated out and (ii) corresponding marginal distribu- 
tions P(9) of the phase 9 when r part of the full dis- 
tribution P(r, 9) is integrated out for the three cases of 
disorder regimes considered above are the (a) weak , 
(b) medium and (c) strong. 

Fig.|| clearly indicate that the phase (9) distributions 
are very different for the three different regimes of dis- 
order parameters considered here. However, the reflec- 
tion coefficient (r) distribution does not show much vari- 
ation with the strength of the disorder parameter. In 
fact, they are quite similar. We will emphasize this point 
again when we later discuss the one-parameter scaling 
theory. The behavior of the probability distribution in- 
dicates that the distribution of the phase has little effect 
on the distribution of the reflection coefficient P(r), and 
hence on the average transport properties. 



C. Distribution of the phase (9) with the strength of 
disorder parameter 2fc£ 

We now study the phase distribution P(9) with the 
strength of the disorder parameter (2fc£). We consider 
here mainly two cases of fixed lengths: (1) I = 1, and (2) 
the asymptotic limit when / — > oo. 



disorder limit P(9) is insensitive to the strength of the 
disorder parameter. Phase distribution is uniform for the 
weak disorder case. It has a perfectly symmetric peak 
centered at 9 = tt for the strong disorder case. In the 
case of intermediate disorder, the distributions are asym- 
metric with respect to the 9 = tt point. For the medium 
disorder strength parameter, P(9) interpolates continu- 
ously between the strong and the weak disorder limiting 
phase distributions. It is clear from the P{9) distribu- 
tions, that the random phase approximation is valid for 
the condition £ > A = 27r/fc, where A is the wave length 
of the incoming wave and £ is the localization length. 
(This has been checked systematically in our numerical 
calculation.) 



2. The phase (9) distribution P(ff) in the asymptotic limit of 
large lengths 

Solution of the FP equation (Eq.|l^) gives the full prob- 
ability distribution P(r, 9) in the (r, #)-space. For larger 
lengths, one can make the approximation: r 1 and 
P{r) w 6(r — 1). In this asymptotic limit, marginal dis- 
tribution for phase P(9) becomes: 

P(d)= [ P(r,9)S{r -l)dr = P(l,9). (26) 
Jo 

Now from Eq.n3 and EqMwe get: 



dP{9) 
dl 



d_ 

09 

2k£ 



d 

(1 + cos 6») — (1 + cos( 
a9 

8P{9) 
89 ' 



P{9) 



(27) 



Fig.^ shows the plot of P{9) in the asymptotic limit for 
different disorder strength parameters (2£;£), as the pa- 
rameter varies from the weak to the strong disorder lim- 
its. Pictures show that in the interval from to 2tt a 
symmetric double peaked distribution around tt in the 
strong disorder regime; a uniform phase distribution in 
the weak disorder regime; and an asymmetric phase dis- 
tribution in the intermediate disorder regime. These 
distributions indicate that the essential features of P(9) 
have not changed with the r ~ 1 approximation, relative 
to the I = 1 case which is shown in Fig.|[ 

D. Discussions on the phase distribution 



1. The phase (9) distribution P(9) for different disorder 
parameters (2fc£) for fixed length 1 = 1. 

Fig.^ shows the plot of P(9) for the fixed length I = 1 
, with different values of the disorder parameters (2fc£). 
Figures show that in the strong as well as in the weak 



We saw that the random phase approximation (RPA) 
is valid for £ > A. Physical meaning of the RPA is that 
the incoming wave has to undergo multiple reflections 
before escaping a localization length and in the process 
the wave randomizes its phase. For the weak disorder 
case the localization length is large. Hence the phase of 
the wave gets randomized. In the other extreme case of 
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the strong disorder, system tries to behave as a perfect 
reflector. Hence, phase of the reflected wave tries to peak 
at 7r (opposite phase with respect to the incoming wave). 
In the regime of intermediate disorder, P{9) distribution 
is disorder specific and has some bias points to peak. 

At this point let us discuss the results of other work- 
ers. Phase distribution for ID systems have been studied 
earlier by Sulem p3| ], Stone et. al . [p4[ , Jayannavar [ p5[ , 
Heinrichs [p6| and Manna et. al. [^7[. Sulem missed the 
random phase distribution in the limit of weak disorder, 
and his calculated phase distributions show peak near 
±7r, and it does not look like a symmetric distribution. 
The study of Stone et. al. showed for the ID Anderson 
model that (i) for the case of weak disorder the phase 
distribution is uniform, (ii) non-uniform for the strong 
disorder, (iii) and pinning of phase distribution near 2tt 
for very strong disorder and large lengths. Their works 
also confirm that the distribution of the phase is insen- 
sitive to the disorder strength in the limit of weak as 
well as strong disorder limits. Manna et. al. show uni- 
form phase distribution from to 2ir for the weak dis- 
order, peaking of the phase distribution near tt for the 
strong disorder. Jayannavar's calculation for the asymp- 
totic (I — > oo) phase distribution shows a uniform distri- 
bution of phase for weak disorder. However, the phase 
distribution for strong disorder shows several peaks that 
do not agree with our results, though we have numeri- 
cally solved the same FP equation. 



E. Marginal probability distribution P(r) of the 
reflection coefficient r with respect to the disorder 
strength parameter 2fc£ 

Fig.|| shows the probability distribution of the re- 
flection coefficient with the disorder parameter strength 
(2fc£) for the sample length 1 = 1. It shows very little 
dependence on the strength of disorder parameter. But 
the distribution certainly has a small spread. 

We can see that though the phase (9) distributions 
P{9) is quite different for different strengths of the dis- 
order parameter (2fc£), the reflection coefficient (r) dis- 
tribution P(r) does not change appreciably with the dis- 
order strength for a fixed length of the sample. 



F. Probability distributions for: g, p, ln(g) and ln(p) 

Once the marginal distribution P(r) is known, the 
probability distribution of any quantity which is a func- 
tion of r, can be easily calculated through the Jacobian of 
the probability transformation. Thus, from the Landauer 
four-probe resistance formula Eq.(|l4|) one can easily cal- 
culate the probability distributions for the resistance(p), 
conductance(g), Inp and Ing , by multiplying with the 
proper Jacobian in to P(r). We can write: 



P(p) = P(r)(l-r) 2 , 
P(g) = P(r)r 2 , 
P(ln(pj) = P(r)r(l -r), 
P(-ln(g)) = P(r)r(l - r). 



(28) 
(29) 
(30) 
(31) 




(Here we have used the same symbol P for before and 
after the transformation.) 

^_shows the plot of P(p ), and 

shows the plot of P(g), for the a fixed length 
10, and for the three different regimes of disorder. 
These distributions show a log-normal form for all disor- 
der strengths. 

Fig.[ll] shows plots of P{ln(p)){= P(-ln(g))) vs ln(p) 
(= —ln(g)), for the fixed length I = 10. These distri- 
butions are Gaussian. (Probability distribution for the 
ln{p) or —ln(g) being a Gaussian implies that the p and 
g obey a log-normal distribution.) 

Fig.[l2| shows the ln(p)(= —ln(g)) distribution for the 
medium disorder with disorder parameter 2fc£ = 1, for 
different lengths /. The distribution matches with the 
Gaussian form quite well even for the case of I = 1. 

We can conclude that the probability distributions for 
p and g are log-normal even for a short sample length (~ 
localization length). They show little dependence on the 
strength of the disorder. The nature of the distribution 
is not affected by the phase fluctuations. 



G. Mean and rms fluctuation of r, g and p 

Figs.|l3|(a),(b),(c), show the plots of the averages and 
the root-mean-squared fluctuations of r, p and g against 
the sample length I for different disorder strength param- 
eters. From the figures one can see < r > first increases 
rapidly and then saturates slowly to r = 1 for the large 
lengths. The average resistance (< p >) shows exponen- 
tial increase and the average conductance (< g >) shows 
an exponential decrease with the sample length I. The 
fluctuations in the case of p and g are always more than 
the average as shown in the figure. However, in case of r 
the fluctuations are always less than the average and are 
finite. At this point we would like to mention that the 
average quantities like average resistance do not admit 
proper treatment by this type of numerical calculation for 
large length. This is because the Landauer four-probe re- 
sistance formula is p — y^— , that is, p is singular at r=l. 
The P(r) never saturates (i.e., no steady state solution), 
but is localized near r = 1 for large lengths. For larger 
lengths, contribution to < p > is dominated by the re- 
values very near to r = 1 (= p — > oo) singularity, and 
numerically it is very difficult to include this long-tail 
p — distribution for the resistance. However, numerical 
calculation of the probability distribution of a quantity 
which is a non-singular function of r, does not pose much 
problem. 
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H. One-parameter scaling theory and the phase 

Here we are re-examining the effect of disorder strength 
and phase distribution on the one-parameter scaling the- 
ory of localization. The scaling theory has been discussed 
in Ref . JH,|2^] . According to the ansatz of the scaling the- 
ory, the dimension-less conductance g (in units of e 2 /irh) 
obeys a universal scaling relation, 



dln(g) 
dln(l) 



(32) 



The above scaling function /3(g) depends only on the 
dimensionality. The scaling relation suggest that the 
fractional change of the conductance per unit fractional 
change of the length depends only on g at that length; 
that is, once g is known at any length scale, one can de- 
rive the conductance for all lengths. The central idea of 
the scaling theory lies in the definition of the /3-function 
which is assumed to have a smooth one-parameter scaling 
behavior with the system size. Now, it is well established 
that these assumptions are in general not correct because 
the scaling theory neglects fluctuations. For the case of 
pure elastic scattering, the resistance/conductance is not 
a self-averaging quantity. Fluctuations increase with sys- 
tem size for insulating regime, and in the metallic regime 
there is UCF, which has a finite non-zero value. These 
issues have been debated by several authors (3tj [33] . 
The points are the following: 

(1) The precise meaning of the scaling function (3(g) is 
not very clear, 

(2) j is a random quantity which depends on the mi- 
croscopic details of the impurities inside the localization 
length, and only for length scales larger than the localiza- 
tion length g may have some universal properties. Only 
defining g, independent of the microscopic details, will 
not specify the average behavior of the sample. 

The question then is how to define a scaling function 
which can accommodate fluctuation effects too. It is 
meaningful to define statistical properties like the proba- 
bility distribution of the resistance/conductance in terms 
of which one can calculate the average properties of the 
system. One can then look for the number (one, or more) 
of parameters needed to specify the probability distribu- 

tion @g. 

There have been several studies 110 Jl6 04J35I1 that show 



that for ID in the weak disorder limit (i.e. within RPA) 
, ln(p) has a normal distribution (we also have discussed 
the same in the earlier Section), in the limit of large 
lengths. 

Now, the expression for the Gaussian distribution for 
logarithm of the resistance will be: 



P(lnp) = 



1 



V2^ 



zexp 



(Inp— < Inp >) 



2a 2 



(33) 



which requires two parameters: the mean < Inp > and 
the standard deviation a. 



It was found |36|] that in the random phase approxi- 
mation (RPA) the nth cumulant C n for the large lengths 
I, 



C n (\n.p) -> la n , 



(34) 



where a n is some constant. It has also been seen by an 
exact solution for the weak disorder case flGl that. 



C\ = a\l and C'2 = a^l 



Eq. (|33|) can now be rewritten as: 



PQnp) 



1 



exp 



(Inp — aiiy 
2ml 



(35) 



(36) 



The above Gaussian distribution has clearly two param- 
eters (al and a2) but can be reduced to a single param- 
eter if the average and the rms fluctuations (variance) 
are universally related. Knowing the probability distri- 
butions P(r) for different lengths and disorder strengths, 
we can check the validity of the one-parameter scaling 
theory for the probability distributions P(lnp). 

In Fig.|l4| we have plotted var Inp vs < Inp > for the 
cases of(l) weak, (2) medium, and (3) strong disorder. 
To obtain these plots we have calculated < Inp > and 
var (Inp), for a chosen value of the disorder parameter 
2fc£, as function of the sample length. Different plots 
thus correspond to the different value of the disorder pa- 
rameter strengths. A One-parameter scaling would have 
these plots coincide. The figure shows that the var (Inp) 
vs < Inp > plots follow approximately the same graph 
for the different disorder strengths and nearly coincide for 
large lengths. To this extent we have a deviation from 
the one-parameter scaling theory in ID. 



VIII. DISCUSSION AND CONCLUSIONS 

We have solved here the ID transport problem for the 
case of Gaussian white-noise disorder numerically. This 
is a nearly complete solution of transport properties of 
the ID Gaussian white-noise random potential that goes 
beyond the conventional random phase approximation 
(RPA) valid only for weak disorder. We have evolved the 
full probability distribution in the reflection coefficient (r) 
and the associated phase (9) space (i.e. (r, #)-space) of 
the complex reflection amplitude R = \/~(r)e ie for a ID 
disordered sample, with different lengths and with dif- 
ferent disorder strengths. For our numerical solution, we 
have taken a fixed initial reflection coefficient ro = .01 for 
all realization of disorder as the Fokker-Planck equation 
for P(r, 9) is singular at r — 0. Gross statistical proper- 
ties of the system are not expected to change with this 
weak extra scatterer. It may, however, affect the very 
sensitive details corresponding to the limit r — * 0, I — ► 0. 
Our numerical work is a systematic study to observe the 
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contribution of the phase fluctuations to different aver- 
ages. 

On the basis of the results obtained, our conclusions are 
the following: 

(A) P(r, 9) and P{6) distributions. 

1. The probability distribution P(r,9) obeys two- 
parameter scaling to include the phase fluctuations and 
merge to one-parameter scaling in the limit of weak dis- 
order. Full distribution of P(r, 9) can be described by 
the length I and disorder parameter 2kt;. 

2. Phase distributions P(9) also have the same two- 
parameter scaling dependence for arbitrary disorder 
strengths; and uniform phase distribution in the limit 
of weak disorder, i.e., one-parameter scaling. 

(B) Phase 9 distribution P{6) in different 
regimes of disorder. 

1. The random phase approximation (RPA) implying 
uniform phase distribution over 2-7T is valid for the condi- 
tion £/A > 1 (£ is localization length and A is incoming 
wavelength), that is, in the weak disorder limit. Phys- 
ically, this means that the wave has to undergo multi- 
ple reflections before it moves through one localization 
length. P(9) is independent of the disorder in the weak 
disorder limit (2/c£ >> 1). 

2. In the strong disorder regime the distribution of the 
phase is perfectly symmetric in the interval from and 
2tt, centering at 9 = n. The distribution is independent 
of the disorder in the strong disorder limit (2kt; << 1). 

3. In the intermediate disorder regime, phase distribu- 
tion is asymmetric about 9 = ir point in the interval 
and 2ir . Also, the distributions are strongly disorder 
dependent. 

(C) Reflection Coefficient (r) distributions P(r): 
The probability distribution for the reflection coefficient 
peaks at r = 1 for large lengths I = L/t; >> 1 
Though probability distribution P{9) for the phase varies 
qualitatively with the variation of the disorder strength, 
the probability distribution for the reflection coefficient 
P{r) does not change qualitatively with the disorder 
strength for a fixed length of the sample. 

(D) Distribution of p, g, ln(p), ln{g): 

For large lengths of the sample (larger than the localiza- 
tion length) the probability distribution for pj g obeys 
the log-normal distribution (i.e., ln(p)/ln(g) obeys the 
normal distribution) for all disorder strengths. An asym- 
metric phase distribution does not affect the nature of the 
probability distribution of the resistance, or the conduc- 
tance, which stays essentially log-normal. 

(E) The mean and the rms fluctuations of p and 

9- 

1. The average resistance increases exponentially with the 
length of the sample. 

2. The average conductance decreases exponentially with 
the sample length. 

3. The rms fluctuations are more than the average for 



both, the resistance and the conductance. 

(F) The One-Parameter Scaling theory and the 
phase 

The one-parameter scaling theory is approximately 
valid for the resistance and the conductance with 
ln{p) or ln(g) as the correct scaling variable for large 
sample lengths. 
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FIG. 8. Reflection coefficient r distribution P(r) against 
the disorder parameter 2fc£ for a fixed sample length 1 = 1. 



FIG. 9. The resistance p distribution P(p) against the 
disorder strength parameter 2k£ for a fixed sample length 
I = 10. 



FIG. 10. The conductance g distribution P(g) against 
the disorder strength parameter 2fc£ for a fixed sample length 
I = 10. 



FIG. 11. Logarithm of the resistance/conductance 
Inp/lng distribution P(lnp)/P(lng) against the disorder 
strength parameter 2fc£ for a fixed sample length I = 10. 



FIG. 12. Logarithm of the resistance/conductance 
Inp/lng distribution P(lnp)/P(lng) against the sample length 
I in the medium disorder regime for a fixed disorder strength 
parameter 2fc£ = 1. 



FIG. 1. Initial probability distribution, i.e., when the sam- 
ple length I — 0: (a) Probability distribution P(r,6). (b) 
Marginal probability distribution P(r). (c) Marginal proba- 
bility distribution P(9). 



FIG. 2. Evolution of the probability P(r, 9) with the sam- 
ple length I in the weak disorder regime for a fixed disorder 
strength parameter 2fc£ = 100 . Plots are for sample lengths: 
(a) I = 1, (b) I = 2, (c) I = 5, and (d) I = 10. 



FIG. 13. Plot of the average and the rms fluctuation ver- 
sus length I for the different disorder strength parameters 2fc£ 
for (a) Reflection coefficient r, (b) Resistance p and (c) Con- 
ductance g. 

FIG. 14. Plot of the var Inp vs < Inp > for different 
disorder strengths 2fc£. 



FIG. 3. Evolution of the probability P(r, 9) with the sam- 
ple length I in the medium disorder regime for a fixed disorder 
strength parameter 2fc£ = 1 . Plots are for sample lengths: 
(a) I = 1, (b) I = 2, (c) I = 5, and (d) I = 10. 



FIG. 4. Evolution of the probability P(r, 9) with the sam- 
ple length I in the strong disorder regime for a fixed disorder 
strength parameter 2fc£ = .001. Plots are for sample lengths: 
(a) I = 1, (b) I = 2, (c) I = 5, and (d) I = 10. 



FIG. 5. Marginal probability distribution P(r) and P(6) 
separately against the sample length I. Three plots are for 
the fixed disorder strength parameters: (a) Weak disorder, 
2fc£ = 100, (b) Medium disorder, 2k£ = 1, and (c) Strong 
disorder, 2fc£ = .001 (corresponding to Figs. 2, 3 and 4). 



FIG. 6. Phase 9 distribution P(6) against the disorder 
parameter 2fc£ for a fixed sample length 1 = 1. 



FIG. 7. Phase 9 distribution P(9) against the disorder 
parameter 2k£ for the asymptotic limit of large length I — > oo. 
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